Emergence of Quantum-Classical Dynamics in an Open Quantum Environment 
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The conditions under which an open quantum mechanical system may be described by mixed 
quantum-classical dynamics are investigated. Decoherence is studied using influence functional 
methods in a model composite quantum system comprising two coupled systems, A and C, interact- 
ing with a harmonic bath with Ohmic and super-Ohmic spectral densities. Subsystem A is directly 
coupled to subsystem C, while C is coupled directly to the bath. Calculations are presented for a 
model where subsystem A is taken to be a two-level system which is bilinearly coupled to a single 
harmonic oscillator C subsystem. The loss of quantum coherence in each subsystem is discussed in 
the extreme non-adiabatic regime where the dynamics of subsystem A is essentially frozen. Sub- 
system C is shown to lose its coherence rapidly, while subsystem A maintains coherence for longer 
time periods since C modulates the influence of the bath on A. Thus, one may identify situations 
where the coupled AC system evolution effectively obeys mixed quantum-classical dynamics. 



I. INTRODUCTION 

Since it is often impossible to isolate a quantum sys- 
tem completely from its environment, the study of open 
quantum dynamical systems is necessary in many circum- 
stances. Techniques used to treat open quantum systems 
have been developed extensively and used in numerous 
applications.Ertl Projection operator techniquescm and 
influence functional methodsa have been actively used 
in chemical physics in studies of electron, proton, and 
exciton transfer processes in the condensed phase and in 
biological systems.LrEj 

In this paper, we investigate some of the circumstances 
under which an open quantum system can be described 
using mixed quantum-classical dynamics. Quantum- 
classical dynamics is_.used to study condensed phase 
many-body systems JlirEj especially in the .-context of 
methods for treating non-adiabatic dynamics JI3c9 In the 
studies presented here, we consider a composite quantum 
mechanical system AC comprising two coupled subsys- 
tems A and C; subsystem C is assumed to be in direct 
contact with a thermal quantum mechanical bath B. Sys- 
tems with this structure can arise in condensed phase dy- 
namics where certain quantum degrees of freedom {A), 
for example, those associated with protons or electrons, 
may interact directly with neighboring solvent molecules 
(C), which in turn interact with the rest of the solvent 
(£?). It is interesting to determine if the composite sys- 
tem may be treated as mixed quantum-classical system 
where the dynamics of subsystem C, which is in direct 
contact with the heat bath, is classical in character while 
subsystem A retains its quantum nature. Some aspects 
of the dynamics of such quantum-classical composite sys- 
tems have been investigated. til 

In order to investigate this problem we study a sim- 
ple model system where subsystem A depends on spin 
degrees of freedom and subsystem C is a single har- 
monic oscillator bilinearly coupled both to subsystem A 



, studied extensively by influence 
2 " The character of its dynamics is 



and the bath. The bath is a collection of independent 
harmonic oscillators. While this is a highly simplified 
model of the realistic systems discussed above, it does 
capture some essential features of real coupled systems 
and is amenable to detailed analysis. Due to its inter- 
action with the bath, the dynamics of subsystem C is 
dissipative and executes brownian motion. The brown- 
ian motion of a quantum particle governed by different 
potential functions and immersed in a thermal harmonic 
oscillator bath has 
functional methods, 
determined by the system-bath coupling and the spectral 
density and temperature of the bath. In the composite 
system we study, the dynamics of the C subsystem os- 
cillator is also influenced by the quantum dynamics of 
subsystem A. The dynamics of subsystem A is more 
complicated. It is also dissipative but its energy must 
be transmitted through subsystem C to the bath. Our 
results on the applicability of mixed quantum-cla.ssical 
dynamics are based on the nature of decoherencecScS in 
the coupled system: when one subsystem behaves quan- 
tum mechanically and the other classically, there must 
be a mechanism making the former decohere slowly and 
the latter quickly. 

The outline of this paper is as follows: In Sec. || we 
specify the model in detail and outline the application 
of the influence functional formalism to it. Although we 
use influence functional methods, similar results can be 
obtained by other methods. Section III considers the 



equilibration of the CB system in the absence of A and 
establishes that this composite system can be described 
by an effective spectral density. The emergence of mixed 
quantum-classical dynamics is investigated in Sec. IV. In 
this section we use the results developed in the previous 
sections, specialized to a two-level system in the strongly 
non-adiabatic regime where the populations change very 
slowly. We study the decoherence of the off-diagonal el- 
ements of the subsystem A reduced density matrix in 
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this regime and show under what conditions one may 
observe quantum-classical dynamics. The conclusions of 
the study are presented in Sec. [v|. The Appendices con- 
tain additional details of the calculations for Ohmic and 
super-Ohmic spectral densities. 



II. GENERAL FORMULATION 

A. Composite system in a thermal bath 

The composite ACB quantum system consists of three 
coupled subsystems: subsystem A is coupled directly to 
subsystem C, while subsystem C is in direct contact with 
a thermal bath B. The ACB system hamiltonian is, 



H = H A + V AC +H C + V CB + A^c + Hi 



(1) 



The thermal bath is taken to be a set of N indepen- 
dent harmonic oscillators with frequencies u> n , masses 
m n , and coordinates and conjugate momenta (q,p) = 
(<ji, qN,Pi, ...,pn) so that the hamiltonian is 



Pn ,1 2 -2 \ 



(2) 



Subsystem C consists of a single harmonic oscillator with 
renormalizcd mass M and frequency u) r with coordinate 
and momentum operators (x,p), 



rr f M 2-2 



(3) 



The interaction potential between C and the bath is as- 
sumed to be bilinear, 



V C 



N 

u = x ^ c n q n , 

71=1 



(4) 



where c n is the coupling constant to the nth oscillator. 
The coupling constants are related to the spectral density 
of the bath by, 



Jb(oj) = 7T - — - — S(lu - U}„). 

* — ' Zra n uj n 



(5) 



We suppose the spectral density has the form 

J fl (w) =r,L0 v e- UJ / h , (6) 

and consider both Ohmic (y — 1) and super-Ohmic 
{v = 3) cases. A counter term AHc has been added 
to the hamiltonian. It depends on c„, m n , uj n , p and x 
and is given by 



r)Ax 2 /ir (v = 1) 

rjAp 2 /M 2 ir + ?yA 3 i 2 /7r (v = 3). 



(J) 



This term is introduced to cancel the shift of the mass 
and frequency of the C subsystem oscillator due to the 
interaction with the bath which will become divergent 
when the frequency cutoff A goes to infinity. As is cus- 
tomary, we consider the renormalizcd quantities after in- 
cluding a counter term as physical observables with spec- 
ified values in our analysis E3 The hamiltonian of sub- 
system A, Ha((j), depends on spin variables a and we 
assume subsystem A is bilinearly coupled to subsystem 
C. Vac = Aai. 

The dynamical properties of interest can be computed 
from the density matrix of the system at time t. In 
the coordinate representation {x, q} of the C subsys- 
tem and bath B, the density matrix takes the form, 
p(x,x' ,q,q' ,t) = (x q\p(t)\x' q'), which is still an opera- 
tor in the A subsystem degrees of freedom. The density 
matrix may be written more explicitly in integral form 
in terms of the propagator 



K(x, q- 1 | x , 9o ; 0) = (x q\e- i6t / h \x 



as 



p(x, x', q, q', t) = dx dx' Q dq dq' Q K(x 1 q; t | x ,q \ 0) 



(8) 



x p(x o ,x' ,q Q ,q' o ,0)k*(x',q';t | x' o ,q' o ;0) . 

(9) 

We are primarily interested in the dynamics of the 
composite AC subsystem under the influence of the ther- 
mal bath. In such a circumstance the reduced density 
matrix p r , obtained by integrating over the bath degrees 
of freedom, is the relevant quantity. Such a reduction 
is especially appropriate if the characteristic time scale 
for the bath is much shorter than those for the A and 
C subsystems. We assume a factorized initial condition 
between the AC subsystem and the bath, 

p(x Q ,x' o ,q o ,q' o ,0) = pac(xq,x' ,0) ® p B (q Q ,q' ,0) , (10) 

and the bath is always taken to be initially in thermal 
equilibrium. Under these conditions the integral form of 
the reduced density matrix at time t is 



p r {x, x',t) = J dx dx' Q J r (x, x'; t \ x Q , x' ; 0)p A c{xo, x' , 0), 

(11) 

where the time evolution kernel J r is given by 
J r (x,x';t | x ,a; o ;0) = / dqdq dq' Q K(x, q;t | x 0l qo]0) 



x PB{qo,q' Q ,0)k*{x',q;t | x' Q ,q' Q ;0). 

(12) 

Given the initial conditions discussed above, this exact 
expression for the reduced density matrix specifies a non- 
Markovian time evolution since the solution at time t de- 
pends on its past history. Approximate Markovian evo- 
lution equations may miss essential features of the quan- 
tum/classical correspondence and tend to underestimate 
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the loss of quantum coherence. In addition, a Marko- 
vian approximation is not generally valid for a harmonic 
oscillator model, except for systems with Ohmic type dis- 
sipation in the high temperature regime. Below we use 
influence functional methods where the exact solution is 
expressed as a path integral in order to deal with the 
nonlocal time evolution. 



B. Influence functional method 

The formulation of the reduced density matrix in terms 
of an influence functional for the composite ACB system 
parallels that -fox. spin-boson models often discussed in 
the literature HEJ Consequently the presentation of the 
generalization to our system will be brief. In order to 
write the reduced density matrix in terms of an influence 
functional, we first introduce a basis of spin functions 
with labels a and a' to represent the states of the A sub- 
system. In this basis we may write eq. (DTI) as, 



p r (cr, a', x, x' , t) — J d<Joda' J dxodx' 

x J r (a, a', x, x'; t | cr , cr'o, ^o, x' Q ; 0)p r (cr o , a' , x , x' , 0). 

(13) 

For the factorized initial condition ( |l0| ) the kernel J r can 
be written in the path integral form as 

J r {a,a\x,x';t | a o ,a' o ,x Ol x' o ;0) = 

p(<t<t') 

/ VoVo' K[o\K*[o'\J a c a (x,x';t | x ,x' ;0), (14) 

where K[a\ is a probability amplitude for subsystem A 
following the path a in the absence of the subsystem C 
and the environment. The evolution kernel Jq 7 for sub- 
system C in the presence of external sources (cr, cr') is 

J^'(x,x';t\x ,x' ;0)= t > VxVx 1 e lS[a ' a ' ' x ' x ' ]/n , 

J(x x' Q ) 

(15) 

The action <S[cr, cr', x, x'\ consists of several contributions 
and can be decomposed as follows: 

S[a, a' ,x, x'] — Sc[x, x'] + Sac[&, x'} 

+AS c [x,x']+Si F [x,x']. (16) 

The actions for subsystem C, plus its counter action, and 
the interaction between subsystems A and C are 

(S C + AS c )[x,x'] = j'ds^x^s) - 

-^'\s) + ^ x%)), 



SacW, a ' i x i x '\ = ^ / ds (x(s)a(s) — x' (s)a' (s)) 
Jo 



where, for notational convenience, we used the bare mass 
Mq and bare frequency ujq as Mq = M and MoUq = 
MljI + 2ryA/vr for v = 1 while M = M + 2tjA/t: and 
Afo^o = MuJ r + 2r/A 3 /(37r) for v = 3. 

The influence action Sif[x,x'] accounts for the effect 
of the bath on C and is given by 

Sif[x, x'] — i ds ds'[x(s) — x'(s)] x 
Jo Jo 

[a(s - s')x(s') -a*(s- s')x'(s% (18) 

where a(t) is a complex kernel whose real a R (t) and 
imaginary a 1 (t) parts, respectively, are given by 

a R {t) = — / dwJs (a;) coth cosut, (19) 

n Jo 2 
1 f°° 

a I (t)=-- duo J b (lo) sin uot. (20) 



We define new variables as R = (x + x')/2, r = x — x' , 
a± = <j±a' ', and write the actions in these new variables 
as 



A f* 

SAc[<r±,R,r] = - J 



ds (cr + (s )r(s)+2a^(s)R(s)). 



{S c + AS C )[R, r] = / ds{M R(s)r(s) - M uj%R(s)r(s)} , 
Jo 

S IF [R,r]=i [ ds [ ds'r(s)a R (s - s')r{s') 
Jo Jo 

-2 [ ds f ds'r(s)a 7 (a- s')R(s'). (21) 



C. Euler-Lagrange equations and solution 



^From eqs. ( f2l|) the Euler-Lagrange equations for R 
and r are 



M R c {s) + M Q ulR c {s) + 2 dsV(s - s')R c (s') 

Jo 

= iU<7 + (s), (22) 



M f c (s) + M u)lr c {s) -2 J ds'a^s - s')r c (s') 

= Act_(s). (23) 



(17) 



The initial and final conditions for eq. ( [22| ) (eq. (|23 
are i?o and R t (r and r t ). If we let the two inde- 
pendent solutions of the homogeneous part of eq. ( p^ ) 
(eq. pg))) be Ui(s)(vi(s)), i — 1,2, with boundary con- 
ditions ui(0) = = 0, u 2 (0) = 0,u 2 (t) = 1 
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(ui(0) = l,vt(i) = 0, u 2 (0) = 0,v 2 (t) = 1), the solu- 
tions of these uncoupled equations can be written as 



R c (s) = RqUi(s) + R t u 2 (s) + -a + (g(s)), 
r c (s) = r Q vi(s) + r t v 2 (s) + /i<7_(g(s)), 

where we have used the notation, 



(24) 



ds'g^(s — s )cr+(s ) 



u 2 (s) ds'g+(t - s')a+(s'), 



a_(g(s)) = / ds'g_(s~s')<j_(s') 



v 2 (s) / ds'g-(t- s')a-(s'). 



(25) 



and defined [i = X/M. The solutions v\ and «2 sat isfy 
the homogeneous part of the backward time equation ( |23| ) 
and are related to u\ and u 2 by i>i(s) = u 2 {t — s) and 
U2(s) = — s). The function g + (s) (g_(s)) also sat- 
isfies the homogeneous part of eq. ( p2| ) (eq. (p3|)) with 
boundary conditions <7±(0) = 0, #±(0) = 1. 

The solutions for gj- are given in Appendix A for 
Ohmic and super-Ohmic spectral densities. From these 
solutions ui. 2 and v\^ 2 can be determined. 



D. Reduced density matrix solution 

Since the potentials in our model are harmonic, an ex- 
act evaluation of the path integral can be carried out. It 
is dominated by the classical solution given in eq. (24). 
The a dependence in eq. (|2~i| ) arises from the back reac- 
tion of the A subsystem on the C subsystem. To separate 
out the contribution from this back reaction, we expand 
S in powers of a and write, 



S[a ± ,R c , r c ] = [R c , r c ] + S« [a ± ,R c , r c ] + S (2) [*±\. 

(26) 

The zeroth order term in a takes the form 
SW[R c) r c ] = (Mux^Ro + Mu 2 {t)R t y t 
- (Mui(0)Ro + Mu 2 (0)Rt)r 
+ i(a n (t)rl + {ai 2 (t) + a 21 (t))r r t 
+a 22 {t)ri) , (27) 



and is ap .influence action for a damped harmonic oscil- 
latorMfel Here 

akl{t) = \l ds f ds'v k (s)a R (s- s')vi(s'), (28) 
* Jo Jo 



for (k,l = 1,2) contains the effects of fluctuations due to 
the bath on subsystem C. The term linear in a is 

S w [a±,R c ,r c ] = A f ds ( Ul (s)R + u 2 (s)R t ) tr_(a) 
Jo 

-\-ifXTQ ds ds' {vi(s)a R (s — s')(7-(g(s'))} 



+ ifir t ds ds'{v 2 (s)a R {s~ s')<r_(g(s'))} 
Jo Jo 

+ ^f ds Ms)v+(s) + ^J* dsv 2 (s)a + (s) . (29) 

This term arises from the interaction between subsystem 
A and the bath, modulated by the trajectory of subsys- 
tem C . The quantum back reaction of subsystem A on 
C induces self-coupling in A, which is contained in the 
last term, 

SW[a ± } = ^J*ds a.(s)a + (g(s)) 



ds I ds' <j-(g(s))a R (s- s')cr-{g(s')). 

(30) 



IjJL 



Using the results above, Jg CT in eq. (|15|) can be written 
in the compact form, 

J a c a '{R u n; t | R , r ; 0) = N(t) exp { -£}, (31) 

where C = R T uf+ ir r &r + o^R + trjr + [a±] and 
N(t) is a normalization factor independent of dynami- 
cal variables. In writing this equation we have intro- 
duced the notation and (a)y = a%j, R T = (Ro,Rt) and 
r 7 " = (ro,r t ), where T stands for the transpose. The 
matrix u is defined as 



"11 «12 \ _ M ( -ui(0) u\{t) 
u 2 i u 22 J ~ y -ua(0) u 2 {t) 

while the vectors a u and a v are given by 



cr_(ui) 
*-(«a) 7 ' 



cr(wi) 

<j{v 2 ) 



(32) 



(33) 



In eqs. (|33]) we have used the symbolic notation 



cr_(ttj) 



dsui(s)<j- (s) 



cr(vi) = / dsw 4 (s)[cr + (s)/2 



+ / ds'ia R (s~s')a^{g(s'))/M] , 
Jo 



(34) 



fori = 1,2. 

Having derived the expression for Jg"' in eq. (|3l|), we 
may use it to obtain the partial Wigner transform of the 
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reduced density matrix—The density matrix in the partial 
Wigner representation^ over the C subsystem degrees of 
freedom is defined by 

p rW (a,a',R,P,t) = ^- [dre- iPr / R 

xp r (a,cr',R+r/2,R-r/2,t). (35) 



Wigner transforming eq. (|13|) and using this definition 
we obtain 



r r{<™') 
p rW {a,a',R,P,t) = da da' / VaVa'K [a]K*[a'] 

J J (ana') 



>1R„(U\j K° a '(R,P;t | R Q ,P o ;0) 
xPrW(Vo,v' ,Ro,Po,0) , (36) 
where the kernel Kg 7 ' (R, P; t \ Rq, P ; 0) is defined by 

K^'(R,P;t | R o ,P o ;0) = -L / drdr e -^ Pr - p ^' h 

A'Ktl J 

xJZ°'(R,r;t | JZo.ro; 0) . (37) 

Inserting the expression for J^f in eq. (|3l|), the time 
evolution kernel Kg 7 for the Wigner function is given 
by 

K^'(R t ,P t ;t | i? o ,P o ;0) 

N{t) f drdro j{-r^+Por +C)/h 



2irh 
N w (t) exp 



(38) 



where Nw(t) = N(t)/2%y\a\ and |a| is the determinant 
of a. The vector SXw = X\y — (X) , with 



Xw = 



Rt 

Pt - \<j(v 2 ) 



(39) 



and 

(X) 



(R) 
(P) 



zl ( «n 1 

U21 V l U l U22 



i?0 

P + A<7(«l) 



The matrix Q is given by 
Q = — x 

"31 



an 



ai2M2i - anu 2 2 



ai2«21 - anU 2 2 anU22 ~ 2ai 2 W 2 lW 2 2 + 022^21 



(41) 



In the next section we shall need Qn an d the explicit 
expression for this matrix element is 



hl(t) = ^ / dwCOth f ^ ) J U U-)±U: I). 



(42) 



where A(w,i) is given by 



A(w,t) = ds ds'g + (t — s) coslu(s — s')g + (t — s'). 
Jo Jo 

(43) 

These results provide the tools needed to analyze the 
conditions under which passage to the quantum-classical 
limit may be carried out. In the next sections we present 
explicit results for both the Ohmic and super-Ohmic bath 
spectral densities. 



III. THERM ALIZATION OF QUANTUM 
BROWNIAN MOTION 

Whether a linear or nonlinear system interacting with 
a bath will eventually reach thermal equilibrium or not 
is one of the fundamental issues in statistical mechanics. 
For a linear quantum brownian motion, from a general 
class of initial conditions, the state of C is-, considered to 
reach asymptotically the gaussian formH In this sec- 
tion, first we will show that the C subsystem evolution 
starting from an arbitrary initial condition possesses an 
asymptotic limit. We then study the influence of such 
a thermalized quantum system on A. For our analysis, 
we shall need the results for a system composed of the 
particle C and the thermal bath B in the absence of A. 
In this case eq. ([36]) takes the form, 

p rW {Rt,P t ,t) = J dR a dP a Kc(Rt,Pf,t\R a ,P o ;0) 

xp rW (R ,P ,0) , (44) 



In the absence of subsystem A, from eq. ( p8[) the time 
evolution kernel for the Wigner function is 

K c (Rt,Pf,t | Ro,Po;0) = Nw^expl-SX^Q-^Xw], 

(45) 



(40) where X^ = (J?, P) and 

<*>-(8)-£(iU)(S)- « 

For a harmonic potential, as a result of the Ehrenfest 
theorem, the center phase space coordinate (X) of this 
distribution function follows a classical trajectory and 
will decay because of energy dissipation into the bath. 
Using the exact expressions for the matrix elements of 
u, obtained from the formulas given in Appendix A, one 
may show that (X) vanishes in the long time limit for 
both Ohmic and super-Ohmic baths. As a result, in both 
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cases, for times long compared to the characteristic re- 
laxation time for this decay, the time evolution kernel 
reduces to 

K c (R t ,P t ; t | Rq, P ; 0) -> N w (t) expl-X&Q^Xw], 
= K c (R t ,Pf,t) (47) 

which is independent of Rq and Pq. Consequently, 
eq. ( [44"| ) reduces to 

PrW {R u P u t)=K c (Rt,Pt;t); (48) 

thus, for an arbitrary initial condition, the Wigner trans- 
formed reduced density matrix is uniquely determined 
by the time evolution kernel Kq- Furthermore, the off- 
diagonal elements of Q vanish in this limit and we obtain 
a Gaussian form for the density matrix whose widths are 
uniquely specified and given by ((Ai?) 2 ) = Qn/2 and 
((AP) 2 ) = Q 22 /2. The width ((AP) 2 (oo)) can be ob- 
tained from eq. (E2I) as 



<(Ai?) 2 M> 



irM 2 



(/a.' COlll ( ^-z I Jb(w)A(w, OO), 

(49) 



We show in Appendix A for the Ohmic case and for 
the weak coupling limit of the super-Ohmic case, that 
J B (w)A(w, 00) = M 2 x'c(u), where 



1 



M (uj 2 - uo 2 + ujj"(uj)) 2 + w 2 7 ,2 (w) ' 



(50) 



is the dynamical susceptibility for C. For the spectral 
density defined in eq. (|^), the explicit form of the fre- 
quency dependent dissipation coefficient 7(0;) is given by 



M-y(u)) 



rjuj 



{u = l) 
(i/ = 3) 



(51) 



in the limit of a large cutoff parameter A, where "/"{uj) 
was absorbed in the mass renormalization for the v = 3 
case. Consequently, we may write ((Ai?) 2 (oo)) as 



(A 2 P(oo)) 



n 

7T Jo 



dtu coth 



(3hu 



XcM- (52) 



Equation (J52) is ,-a statement of the fluctuation- 
dissipation relation.EJ The form in eq. (^2|) suggests that 
the combined system CB is now in thermal equilibrium 
with spectrum specified by x'c(<jj) but ^ e system C itself 
is not in equilibrium. The fluctuation-dissipation relation 
is violated at finite t while the system is still far from equi- 
librium and holds only asymptotically. This result gen- 
eralizes the result previous]* obtained by Caldeira and 
Leggett for an Ohmic bath.cj 

When the CB system is in thermal equilibrium, we 
show in Appendix |b| using linear response theory that 
Jcb{^) = A 2 x'c( w ) : the effective spectral density is pro- 
portional to the dynamical susceptibility of the C system. 
If we use the relation 



Jb(cj) = Mw7'(w), 



(53) 



for the bath spectral density assumed in eq. (g), we see 
that the effective spectral density in the combined system 
CB can be written as 



Jcb(oj) 



M 2 {lo 2 -lu 2 + ^"{lu)) 2 + J 2 (lo) ' 



(54) 



This relation can also be obtained by solving, the classical 
equations of motion as suggested earlierErE§E3 owing to 
the linear coupling assumed between the C system and 
the bath. 



IV. EMERGENCE OF QUANTUM-CLASSICAL 
DYNAMICS VIA DECOHERENCE 



The investigation of the emergence of quantum- 
classical dynamics presented here will be restricted to 
cases where the dynamics of subsystem A occurs on time 
scales which are very long compared to those that char- 
acterize subsystem C. We suppose that subsystem A is a 
two- level system with a — ±1 and focus on the extreme 
non-adiabatic regime where the population dynamics of 
subsystem A is essentially frozen on the time scales of in- 
terest so that pa(— 1j — lj t) = /0^(— 1, — 1, 0), pa(1, 1, t) — 
Pa(1, 1)0). We further assume that subsystems A and C 
are decoupled initially and use the factorized initial con- 
dition, 



Pac(xq,x' o ,0) = p A (0) ® pc(x ,x' ,0). 



(55) 



This is a reasonable assumption if they are weakly cou- 
pled. 



A. Decoherence in A interacting with an equilibrium 
CB bath 



If the interaction between the A system and the com- 
bined system CB is turned on after an equilibrium state 
is reached, the effect of CB on A is given simply by the 
harmonic bath directly coupled to the A subsystem but 
with its spectrum characterized by the effective spectral 
density. Under these conditions, the off-diagon al p art 
of the density matrix can be obtained from eq. ( |l3| ) by 
dropping the dependence on (x,x'), omitting the inte- 
grals over spin variables and by making the variable re- 
placements x — > a and Jb — * Jcb in eq. ([l8]) for the 
influence functional. We find, 



p4-l,l,t) = e- D i O) (V(-l,l,0), 



(56) 



where D^\t) is a decoherence factor for A and is defined 



as 



G 



4 



duJcsiu) coth 
cos(wt) 



(57) 



Thus, as mentioned above, the dynamics of the compos- 
ite system is equivalent to that of a two-level spin-boson 
system coupled to a harmonic bath, with an effective 
spectral density. 
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FIG. 1. The temporal behavior of the off-diagonal ele- 
ment of the reduced density matrix pa{— 1> 1> t) f° r subsys- 
tem A is plotted in panel (a) versus time. The initial value 
Pa(— 1) 1,0) = 1. Parameters are: (3 = 0.66, A = 500, ui r — 1, 
7 = 0.3, M = 1 and h = 1. The same quantity is shown in 
panel (b) using a logarithmic time scale. In both figures the 
solid line corresponds to A under the influence of CB with the 
effective Ohmic spectral density, while the dashed line shows 
the influence of Ohmic bath B without C . (1) A 2 = 0.1 and 
(2) A 2 = 0.3. 

Using these results we may examine the coherence in 
A under the influence of the effective bath CB and the 
coherence of A under the influence of the bath B sepa- 



rately. In Figs. 1(a) and (b), the off-diagonal part of the 
density matrix of subsystem A is plotted for an Ohmic 
environment. 

The solid line corresponds to our system while the 
dashed line is the ordinary spin boson model without C. 
The time evolution of pa{— T 1>£) f° r f ne ordinary ordi- 
nary spin boson model is characterized by (1) an initial 
period where subsystem A system has not yet felt the 
existence of the environment and its coherence is main- 
tained; (2) a quantum regime for t > 1/A where the 
system begins to interact with the vacuum fluctuations 
of the environment and, finally, (3) a thermal regime for 
t > Ti(3 where the effects of thermal fluctuations on sub- 
system A have set in. For 1/A < t < h/3, only vacuum 
fluctuations interact with the system. 

By contrast, the existence of the intermediate subsys- 
tem C changes the evolution of Pa(— 1, 1, t) significantly. 
Subsystem C is characterized by its harmonic oscillator 
frequency and mass. Fluctuations originating in the bath 
B are modulated by C through these parameters. For a 
large mass M , the effective modes of the combined CB 
system are concentrated in the neighborhood of uj r and 
the decoherence behavior of subsystem A is governed by 
these modes. Therefore, after the initial period, decoher- 
ence will begin to occur at t ~ l/u> r . 

Although some systems are approximately character- 
ized by an Ohmic spectral density, generic systems will 
have non-Ohmic spectral densities. The density of states 
typically varies as uo u with v depending on the spatial 
dimension D. We consider the 3D super-Ohmic case 
with v = 3. Such- an environment is relevant for the 
study of polaroids ,tl!] macroscopic magnetization tunnel- 
ing in crystalsjlll and radiation damping of atomsEa A 
super-Ohmic environment affects subsystem A evolution 
on short time scales more significantly than an Ohmic 
environment. A sub-Ohmic environment induces non- 
trivial long time behavior. Results for the super-Ohmic 
case with v = 3 are shown in Figs. || (a) and (b). 

For the super-Ohmic high dimensional environment, 
the vacuum fluctuations play a more significant role than 
for an Ohmic environment, reflecting the fact that there 
are larger number of high frequency modes in the higher 
dimensional environment. The same effect is also respon- 
sible for the large difference in the behavior of correlation 
functions of the system in quantum and classical baths 
when the bath is super-OhmicJj The dynamics of A in 
the presence of the intermediate subsystem C is worth 
noting. For the present choice of parameters, pa{— 1; h t) 
asymptotically reaches a non-zero value indicating that 
quantum coherence in the A system will never be lost 
completely. The origin of this behavior can be under- 
stood in terms of the time-dependence of the diffusion 
coefficient (See Ref£3). System C initially executes brow- 
nian motion with a time-dependent diffusion coefficient 
as a result of the non-Markovian time evolution of the 
reduced density matrix. This time-dependent diffusion 
coefficient of C, which also determines the decoherence 
rate, exhibits a rapid increase at early times t > 1/A and 
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then asymptotically vanishes. This initial increase can 
be large enough to wash away the quantum coherence 
of C directly coupled to the bath. Owing to the mod- 
ulation effect from C, these high frequency modes are 
filtered out and do not directly affect A if ui r << A. The 
late time value of the diffusion coefficient is too small 
to eliminate the quantum coherence from subsystem A 
completely Hence, the non-Markonvian nature of the 
density matrix evolution for the super-Ohmic case is re- 
sponsible for the significant difference in the dynamics of 
coherence between the A and C subsystems. 
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FIG. 2. Plot of pa{— 1, versus time (a) and versus 
log(t) (b) for a super-Ohmic environment with v — 3. Pa- 
rameters and labeling are the same as in Fig. The dashed 
line in panel (a) for super-Ohmic bath B without C overlaps 
the vertical axis. The behavior is visible in panel (b) using 
the logarithmic time scale. 



B. Decoherence in A for arbitrary C initial 
conditions 

Thus far we have considered only the asymptotic limit 
where CB has reached equilibrium and the fluctuation- 
dissipation relation holds. It is interesting to consider the 
coherence of A under the influence of intrinsic dynamics 
of C starting from arbitrary C initial conditions. As an 
example of this situation, we assume that C is initially 
in a gaussian coherent state given by the wave function, 



tpxp(x) 



1 



(Tre 2 ) 1 /' 



■ exp 



(x — xY 



px 



(58) 



with width e = l/ A /ZZv where x and p are parameters. 
The C subsystem density matrix at the initial time is 
then given by 



pc(x ,x' ,0) 



1 



: exp 



(i?o - xf r ; 



o 
4e 2 



.pro 



(59) 



We may now use this expression for pc(xo, x' , 0) in the 
factorized initial condition for AC given in eq. ( pq ) . The 
off-diagonal element of reduced density matrix for subsys- 
tem A in the non-adiabatic limit can then be determined 
from eqs. jf^ ) and ( |i~4| ) for p r (a,<r' ,x,x' ,t) for AC by 
integrating over the C variables. The result is, 



Pa(— 1, 1, t) = / dxdx' 



dxodx 

x J c 1:1 (x, x',t\x , x' , Q)pc{x , x' , 0) 

xpA(-l,l,0). (60) 

Changing integration variables {x,x') — > (R,r) 
and (xq,x ) — > (Ra,ra) and using eq. ( |3l| ) for 
J^ 1,1 (i?, r, t\Ro,ro, 0) we may write this as 

p A {-l,l,t) = J dRdr [ J dR dr 

xN(t)exp{^L}pc(Ro,r ,0) 



xp A (-l,l,0) . 



Making use of the expression for C evaluated for a = 
and a 1 = 1 and carrying out the integrations we find 

PA (-1,M) = p A {-l,l^)e- DA{t) 

x exp 



(61) 

= —I 



x exp 



j:(f R (t)x + f P (t)p) 



1 



~(fl(t)/Mtj r + fUt)Mw r ) 



where 

/*(*) 
and 



ds^-. 

o u 21 {s) 



fp(t) = A 



ds — 

o "21 W 



(62) 



(63) 
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D A (t) 



2A 2 



irhM 2 



dio Jb(cj) coth 



f fifiuj 



A A {uj,t). 



(64) 



The explicit expression for A A {uJ,t) can be computed 
from the solutions presented in Appendix A by a lengthy 
but straightforward calculation. We find, 



A A {uj,t) 



1 



(n 2 r - uj 2 ) 2 + 4T 2 uj 2 
2 1 — cos(wi) 



1 — cos(wt) 



nn 2 uj 

x - e~ rt cos(Qi)) - cjr e ~ r * sin(fl*)) 

2 sin(wt) 

x (2rfi(l - e~ rt cos(fii)) + (ft 2 - r 2 )e~ r * sin(m)) 

+ tm r ^ 2 + 4r2 ) o2(1 - co « e_rt ) 2 

-2rfi (uj 2 + 2T 2 - 2ft 2 ) 

x (1 - cos(ftt)e _r *) sin(Ot)e _r * 

+ (^ 2 r 2 + (ft 2 - r 2 ) 2 ) sin 2 (fti) e - 2rt }] . (65) 



Equations (p2|)-(|65|) are our main results. The first term 
in eq. ( |62| ) is the initial condition for A, the second term 
accounts for the decoherence arising from thermal fluctu- 
ations of the bath mediated by C. The third term arises 
from the initial condition for C and the last term is re- 
sponsible for the decoherence of A due to the averaged 
damped oscillatory motion of C . 

The modulus of the off-diagonal reduced density ma- 
trix element for A is shown in Figs. ||(a) and (b) as a 
function of time. The solid line shows the second term 
eq. (|62j), S = er DA ^ l \ which gives the decoherence due 
to the quantum-back-reaction-induced self-interaction of 
subsystem A. The dashed line shows the last term, 
L = ex^[-\(f%(t)/Muj r + f%(t)Muj r )], the decoher- 
ence due to the motion of C. The latter effect is not 
strong enough to eliminate coherence of A completely. 
Since, initially, the CB is far from equilibrium, those 
modes that affect A are not necessarily near resonant 
modes around ui r . Consequently, the decoherence of A 
is more rapid in this case than in the case of evolution 
from the thermal equilibrium initial condition studied in 
Sec. IV. A. This tendency is more evident in a super- 
Ohmic bath than in an Ohmic bath. With generic initial 
conditions for C, A is under the influence of the nonequi- 
librium bath CB, which is no longer equivalent to the 
effective thermal bath with the spectral density Jcb- 
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FIG. 3. Plot of contribution to the modulus of Pa{— 1> 1> t) 
versus time for an Ohmic bath (Fig. |^(a)) and for a super 
Ohmic bath (Fig. ^(b)). Subsystem C is in a coherent state 
initially. In both figures the solid line corresponds to the sec- 
ond term S, while the dashed line corresponds to the last term 
L in the right hand side in eq.(^). A 2 = 0.1 for SI, LI and 
A 2 = 0.3 for S2,L2. Parameters values are the same as in 
Fig. % 

Now let us consider the situation where the decoher- 
ence of C is fast and that of A is slow. This may occur at 
high temperatures and for weak coupling between A and 
C. In such a case, the back reaction from A on C may be 
neglected. In this circumstance, we may consider the de- 
coherence of C in the absence of A. Quantum brownian 
motion and decoherence of a damped harmonic oscillator 
has been studiedi^reviously.ta'l^irEZl Using the result in 
eq. (3.9) of Ref.jHj the reduced density matrix may be 
approximated by, 

Pr (a, a', x, x', t) ~ e-^- x '^ Dc ^p r (a, a', x, x', 0), (66) 
where Dc(t) = / * ds J Q S ds'a R (s') costs' is the decoher- 
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ence factor for C. 
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FIG. 4. Plot of the decoherence factor versus time. The 
solid line corresponds to the logarithm of the modulus of the 
density matrix pa in eq.(^) with the inverse sign, while the 
dashed line is Dc- System parameters are: x — x' = 1, u r = 1, 
7 = 0.3, M = 1, % = 1. Panel (a) is for the Ohmic bath with 
/3 = 0.002, A = 500. Panel (b) is for the v = 3 super-Ohmic 
bath with (3 = 0.02, A = 20. 

In Fig. ||, the decoherence factors for A and C are 
plotted. For an Ohmic bath (Fig. 0(a)), initially, the 
quantum coherence is lost faster C than in A; hence, the 
A subsystem behaves more quantum mechanically than 
the C subsystem during this initial period. Figure |J(a) 
shows that there is a crossover and for longer times, 
t > l/u) r , subsystem A experiences stronger decoherence 
than C. The decoherence factor for C varies linearly with 
t indicating that the dynamics of C asymptotically ap- 
proaches Markovian evolution in the high temperature 
Ohmic bath. For a super-Ohmic bath (Fig. (|(b)), a sub- 
stantial portion of the quantum coherence in C is lost 
during the initial period due to large fluctuations coming 



from the bath, while A retains quantum coherence for a 
long time. The dynamics of the coupled system can be 
approximated by mixed quantum-classical evolution for 
long times after the short transient period. 



V. CONCLUSION 

An appreciation of the conditions under which a 
quantum mechanical system may be approximated as a 
quantum-classical system is essential in many applica- 
tions to the dynamics of many-body systems. Adopting 
an open quantum system point of view has provided a 
natural way to explore this issue based on the decoher- 
ence of the superposition of quantum states into a statis=. 
tical mixture under the influence of the environment .L3'E3 
Indeed, considerations of decoherence have played an im- 
portant role in discussions oiLschemes for the simulation 
of condensed phase systems O 

Using influence functional methods, we have shown for 
our simple model system that the decoherence time scales 
that characterize the A and C subsystems can differ sig- 
nificantly. In particular, in the limit of nonadiabatic dy- 
namics we have identified the following three regimes: (1) 
the full quantum regime where both the A and C sub- 
systems behave quantum mechanically; (2) the quantum- 
classical regime where subsystem A maintains coherence 
owing to its indirect coupling to the bath, while C has 
lost its coherence and behaves effectively classically; (3) 
the classical regime where the quantum coherence of both 
the A and C subsystems has been lost and the compos- 
ite AC subsystem exhibits effectively classical dynamics. 
The different roles that the environment plays in its ef- 
fect on two subsystems is responsible for the separation 
of decoherence time scales. 

We also saw that a different choice of initial conditions 
will modify above picture. When the initial condition 
for C is chosen to be a pure gaussian coherent state, CB 
can no longer be considered to be a thermal bath and the 
coherence of A will be lost in a shorter time scale than 
times given by the inverse characteristic frequency of C. 

Quantum/classical correspondence in nonlinear sys- 
tems, in particular, in chaotic systems, has many non- 
trivial featureslij The quantum open system approach 
adopted-, in this paper can also shed light on this prob- 
lem.^' a While the model system studied here is very 
simple and the extreme non-adiabatic regime is only one 
limiting case to consider, our results nevertheless have 
provided some insight into the emergence of quantum- 
classical dynamics and should be useful in the study of 
the quantum dynamics of more complex systems. 
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APPENDIX A: SOLUTIONS OF EQUATIONS OF 
MOTION 

In this Appendix we give some details of the solutions 
of the equations of motion discussed in Sec. [I C which 
are used in the calculations presented in Sec. IV. We 
require the functions iti,2(s) (and Ui,2(s)) which are the 
solutions of the homogeneous parts of the Euler-Lagrange 
equations (|2^) and (|23|) . These functions can be found if 
ff±( s ); which are also solutions of the homogeneous parts 
eqs. (|22|) and fl23j), are known, since 



M • i s 9+( s ) ■ ,.\ i \ 9+{ s ) , A1 s 

ui(s) = g+(s) —g + (t), u 2 (s) = — j-r, (Al) 



9+{t)- 



9+{t) 



in order to satisfy the boundary conditions on these func- 
tions. Using the spectral density for an Ohmic bath, the 
equation that g±(s) satisfies is 



g±(s)±2Tg ± (s)+u 2 g±{s) =0 



The solution of eq. (A2) is 

, . sin fts 



9±{ s ) 



n 



(A2) 



(A3) 



where T — 7 and ft 2 = co 2 — j 2 . 

Using the spectral density for a super-Ohmic bath with 
v = 3, the equations of motion for g±(s) can be written 
as 



'g±{s)+uj 2 r g±{s)T^l 9± (s) =0 



(A4) 



Thus, the equation of motion has the form of the 
Abraham-Lorentz equation.E3 This equation is an ap- 
proximation to the full equation of motion with non-local 
time dependence. As in the electromagnetic field case, 
only the physically relevant roots of the characteristic 
equation must be retained. These are r ± ift with 



1 [D 1 
r = — — + 1 

67 V 4 D 

where 



67 V 4 D 



(A5) 



D = 2 (l + 54( 7 uv) 2 + VW8jLu r v / 27(juj r ) 2 + lj 



1/3 



Retaining the physically relevant roots of the character- 
istic equation, the solution for g+(s) for the super-Ohmic 
case has the same form as eq. (A3) but with the above 
values of T and ft. In the weak coupling limit that we 



consider in our calculations, T 



and ft 



In our model, eq. (A4) can be reduced to a second or- 
der ordinary differential equation by differentiating the 



homogeneous part of the Euler-Lagrange equations with 
respect to time as 



9± (s) = ~Lu 2 g ± (s) + 0( 7 ) 



(A7) 



which may be substituted into eq. (A4) to give ( |A2| ) with 
r = ~fuj 2 and ft = uj r in the leading order approximation 
as indicated above. Altho ugh one can improve the ap- 
proximation leading to eq. ( A2) for the super-Ohmic case 
to arbitrary higher order in 7, we restrict our study to 
the leading order in this parameter. 

Since the u± t 2 and v±.2 solutions are now known, we 
have all the information needed to compute the various 
quantities necessary to obtain the numerical results. In 
particular, the expression for the matrix elements of a in 
eq. (pj^ipi a general environment can be found. (See 
RefsHESEH for an Ohmic bath.) In our computations, 
we need fflii(i) which has the form, 



a n (t) 



2*g+(t) Jo 



dujJ B {u) coth 1 — —A(w,t) , (A8) 



where A(w,i) was defined in eq. (|43[) . The evaluation of 
A(uj,t) is straightforward and leads to, 



A(w,t) 



1 



[2ft 2 + , 



-rt 



2ft 2 (ft 2 ~uj 2 ) 2 +AT 2 uj 2 
x [2cos(wi){(r 2 +lu + lu_) cos(ftt) - 2rftsin(fti)} 
- (r 2 + uj 2 _ ) cos(uj + t) - (r 2 + uo 2 + ) cos(w_t)] 



-2Ft 



2rftsin(2ftt) 



(r 2 +UJ+OJ-) cos(2fti)}] 



(A9) 



where u>± = to ± ft and ft 2 = ft 2 + T 2 . Using this result, 
in the asymptotic limit we obtain, 



J B (u;)A(u), 00) 



{ 2M^ui 
2M~/u 3 



(l/ = l) 

(i/ = 3) 



In writing this equation, we used the fact that the dis- 
tribution J b (uj)A(ui, 00) is highly peaked around u r for 
small 7. 

Similarly, for long times we may write the elements of 
a as, 



an(f) 



2J2Yt 



ft 2 e 



27r sin fit 



dujJ B (uj) coth A{ui, 00) 



a 12 (t) = a 2 i(t) 



(A6) 



fte rt 
27r sin fti 



( T -ftcot(fts))A(w,oo) 



a 2 2(t) 



1 

2^ 



duj J b (lu) coth 



du>J B (u>) coth 



(3huj 



(3huj 



{ lu 2 +(r- ftcot(fts)) 2 } A(w,oc). 



(A10) 



Using these results we can compute the asymptotic prop- 
erties of Q needed in the calculations presented in the 
text. 
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APPENDIX B: EFFECTIVE SPECTRAL 
DENSITY 

In this Appendix we derive the relation between the 
effective spectral density for the CB system and the dy- 
namical susceptibility of subsystem C using an argument 
that by-passes the actual diagOBsJiaation procedure. A 
similar argument is given in RefsE3'E3 From the Hamilto- 
nian given in eq. (Q), we can write the equation of motion 
for C with a harmonic potential in the form, 



(Bl) 



in the complex Fourier representation with a susceptibil- 
ity function xc(w). The force Fc(u>) — —dVAc{u)/dx = 
— \<j{uj) is the external force from subsystem A acting 
on C . On the other hand, the equation of motion for 
subsystem A gives 



(B2) 



with Fa(uj) = —dVA(cr)/da. Here we assumed Ha has 
the form p A /2mA + Va(ct) for simplicity. Our argument 
does not depend on the form of Ha as clear from the 
context. Combining with eq. (Bl), we have 



(B3) 



Now, suppose we have already diagonalized last three 
terms of eq. (|lj) and replaced them with N+l - harmonic 
oscillators, 



He + Vcb + Hb 



H' b 



N 
n=0 



'2 



2m' 



). (B4) 



We write the interaction term Vac m terms of new coor- 
dinates q' n 



Vac -» V' A 



ab 



N 
n=0 



(B5) 



Assuming that the system described by the Hamiltonian 
in eq. (Q) reaches thermal equilibrium in the asymptotic 
limit, we can replace the effect of the combined system 
CB by the equivalent thermal bath H' B at the same 
temperature. Although this diagonalization procedure is 
straightforward, for our purpose, we only need the form 
of the equation of motion for a expressed by the dissipa- 
tion coefficient 7cb(w) for the effective bath eq. (B4), 



[-m A L0 2 - im A LO^cB{^)\ = F a (uj). 



Comparing the above with eq. (B2) we have 
A 2 xcM = im A uJlcB{u). 



(B6) 



(B7) 



Recall that the spectral density is related to the dissipa- 
tion coefficient by 



Jcb(u) = m A uJ^ CB (u}), (B8) 
for real ui. From this the desired relation follows: 

J cb (lu) = A 2 x'cM- (B9) 
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